A practical solution to the sign problem at finite theta- vacuum angle 



Takahiro Sasaki/ Q Hiroaki Kouno, 2 '0 and Masanobu Yahiro'd 

'Department of Physics, Graduate School of Sciences, Kyushu University, Fukuoka 812-8581, Japan 
2 Department of Physics, Saga University, Saga 840-8502, Japan 
(Dated: August 3, 2012) 

We propose a practical way of circumventing the sign problem in lattice QCD simulations with a theta- 
vacuum term. This method is the reweighting method for the QCD Lagrangian after the Ua (1) transformation. 
In the Lagrangian, the P-odd mass term as a cause of the sign problem is minimized. In order to find out a good 
reference system in the reweighting method, we estimate the average reweighting factor by using the two-flavor 
NJL model and eventually find a good reference system. 
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I. INTRODUCTION 

Phenomena based on strong interaction have shown that 
charge conjugation C, parity P and time reversal T are good 
symmetries of nature. This means that quantum chromody- 
namics (QCD) should respect any combinations of the dis- 
crete symmetries. Among the discrete symmetries, the CP 
symmetry is not necessarily respected in QCD due to the in- 
stanton solution [1,2]. The instanton solution allows the QCD 
Lagrangian Cqcd to have a ^-vacuum term. The resulting 
Lagrangian is described as 
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in Euclidean spacetime, where F" is the field strength of 
gluon. The vacuum angle 9 is a periodic variable with pe- 
riod 2ir. It was known to be an observable parameter yfl. The 
Lagrangian Cqcd is invariant under the combination of the 
P transformation and the parameter transformation 9 —> —9, 
indicating that the P and CP symmetries are preserved only 
at 9 = and ±7r; note that 9 = —ir is identical with 9 = ir. In 
the 9 vacuum, therefore, we must consider the P and CP vio- 
lating interaction parameterized by 9. Theoretically we can 
take any arbitrary value between — ir and ir for 9. Never- 
theless, it has been found from the measured neutron electric 
dipole moment 0] that |0| < 10~ 9 Jl-01- Why is 9 so small in 
zero temperature (T)? This long-standing puzzle is called the 
strong CP problem; see for example Ref. [8[] for the review. 

Around the deconfinement transition at T = T4, there is 
a possibility that P-odd bubbles (metastable states) arise and 
thereby regions of nonzero 9 are generated [9]. Thus 9 can be- 
come a function depending on spacetime coordinates (t, x). If 
P-odd bubbles are really produced at T ss Aqqd, P and CP 
symmetries can be violated locally in high-energy heavy-ion 
collisions or the early universe. This finite value of 9 could be 



a new source of large CP violation in the early universe and 
a crucial missing element for solving the puzzle of baryogen- 
esis. 

In the early stage of heavy-ion collision, the magnetic field 
is formed, and simultaneously the total number of particles 
plus antiparticles with right-handed helicity is deviated from 
that with left-handed helicity by the effective 9(t, x). In this 
situation, particles with right-handed helicity move opposite 
to antiparticles with right-handed helicity, and consequently 
an electromagnetic current is generated along the magnetic 
field. This is the so-called chiral magnetic effect IU0l - [l3ll . 
The chiral magnetic effect may explain the charge separa- 
tions observed in the recent STAR results HI. Hot QCD with 
nonzero 9 is thus quite interesting. 

For zero T and zero quark-number chemical potential (/1), 
some important properties are showed on P symmetry. Vafa 
and Witten proved for 9 = that the vacuum is unique and 
conserves P symmetry lfl5ll . This theorem does not preclude 
the existence of P-odd bubbles. At 9 = tt, P symmetry is 
also preserved as mentioned above, but it is spontaneously 
broken lll6l Il7ll . The spontaneous violation of P symmetry 
is called the Dashen mechanism flrjll . Although the mecha- 
nism is a nonperturbative phenomenon, the first-principle lat- 
tice QCD (LQCD) is not applicable for finite 9 due to the 
sign problem. The mechanism at finite T and/or finite /1 
was then investigated with effective models such as the chiral 
perturbation theory IU8I - I23I1 . the Nambu-Jona-Lasinio (NJL) 
model l24U27tl and the Poly akov-loop extended Nambu-Jona- 
Lasinio (PNJL) model JH-SS. 

In the previous work [30], we proposed a way of minimiz- 
ing the sign problem on LQCD with finite 9. The proposal is 
as follows. For simplicity, we consider two-flavor QCD. The 
QCD Lagrangian (Q~|i is transformed into 
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C QC D=q'M(9)q' + -F; v F« v 
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with 



M (9) = lv D v + m cos (9/2) + mi~f 5 sin (9/2) (3) 



„, — , — , by using the Ua(1) transformation 
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where the quark field q — (q u ,Qd) has been redefined by the 
new one q'. The determinant M. (0) satisfies 



det M(6) = [detM(- 



(5) 



indicating that the sign problem is induced by the P-odd (0- 
odd) term, mi^ sin (0/2). The difficulty of the sign problem 
is minimized in ©, since the P-odd term with the light quark 
mass m is much smaller than the dynamical quark mass of 
order /Iqcd- Actually, it was found that the P-even conden- 
sates a j — (q'fq'f) is much larger than the P-odd condensates 
Vf = (<7/275<7/)- The P-even condensates little change even 
if the 0-odd mass term is neglected. We then proposed the fol- 
lowing reweighting method. The vacuum expectation value of 
operator O is calculated by 



(O) = J VAOdctM(6)e- s ° 
= J VAO'detM Tei {6)e- s « 
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with the gluon part S g of the QCD action and 

O' = R{6)0, 
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where R(9) is the reweighting factor and det A4 Ie f(0) is the 
Fermion determinant of the reference theory that has no sign 
problem. The simplest candidate of the reference theory is the 
theory in which the 0-odd mass is neglected. We refer this 
reference theory to as reference A in this paper. As discussed 
in Ref. lT30h - reference A may be a good reference theory for 
small and intermediate 6, but not for large 9 near tt. In refer- 
ence A, the limit of 9 = tt corresponds to the chiral limit that 
is hard for LQCD simulations to reach. 

The expectation value of R(9) in the reference theory is 
obtained by 



(10) 



where Z (Z rc [ ) is the partition function of the original (refer- 
ence) theory. The average reweighting factor (R(9)) is a good 
index for the reference theory to be good; the reference theory 
is good when (R(9)) = 1. 

In this paper, we estimate (R(9)) with the NJL model in or- 
der to find out a good reference theory. We find that reference 
A is good only for small, so propose a good reference theory 
that satisfies (R(6)) « 1. 

This paper is organized as follows. In Sec.[n] we recapitu- 
late the two-flavor NJL model and show how to calculate the 
pion mass and (R(9)) for the case of finite 9. Numerical re- 
sults are shown in Sec. [HI] Section lTVl is devoted to summary. 



H. MODEL SETTING 

The two-flavor NJL Lagrangian with the 0-dependent 
anomaly term is obtained in Euclidean spacetime by 



£ = q(lud v + m )q - G^y] [{qr a q) 2 + [qi^T a q) 2 ] 

- 8G 2 [e w detq R q L + e^de^R,] (11) 



where mo is the current quark mass satisfying mo = m u — 
rrid and To and r a (a = 1,2, 3) are the 2x2 unit and Pauli 
matrices in the flavor space, respectively. The parameter G\ 
denotes the coupling constant of the scalar and pseudoscalar- 
type four-quark interaction, while G2 stands for that of the 
Kobayashi-Maskawa-'t Hooft determinant interaction I2I I3II1 
where the matrix indices run in the flavor space. 

Under the Ua(1) transformation (@}, the quark-antiquark 
condensates, {&} = {a, r/, flj, 7Tj}, are transformed as 

a = qq = cos((9/2)cr' + sm(0/2)?/, (12) 

r\ = qij 5 q = -sm(6/2)a' + cos(0/2)r/, (13) 

Oi = qnq = cos(0/2)a- + sin(0/2)7r-, (14) 

7Tj = qi^nq = -sin(0/2)a- + cos(0/2)tt,', (15) 

where the condensates {^} = {a 1 , rf , a\, 7r^} are defined by 
the same form as {<f>i} but q is replaced by q' . The Lagrangian 
density is then rewritten with q' as 
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q'M + m{6))q' - d [{q'r a q') 2 + (q'i^T a q'f] 

(16) 



a—0 



- 8G 2 [det<? K qi + det^gJJ 
= ~q\l v d v + m{9))q' - G + [(q'q') 2 + (q 1 l^rq'?] 

- G_ [{q'rq'f + [q'fq'f] , (17) 
where G± = G\ ± G2 and 

m(0) = m o cos(0/2) +m o i7 5 sin(0/2). (18) 

Making the mean field approximation and the path integral 
over the quark field, one can obtain the thermodynamic po- 
tential fl (per volume) for finite T: 



Q = U-AN C 
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with 



E± = \Ip 2 + C±2VD, 

C = M 2 + N 2 + A 2 + P 2 , 

D = (MA + NP) 2 + (Ax Pf > 0, 

M = m o cos(0/2) - 2G+cr', 

N = m o sin(0/2) - 2G-T)', 

A = -2G-a', P = -2G+7?', 

A = Vl-A, p = VP-P, 

U = G + (a' 2 +7f' 2 ) + G^(r 1 ' 2 +a' 2 ), 



(19) 



(20) 
(21) 
(22) 
(23) 
(24) 
(25) 

(26) 
(27) 



3 



where the momentum integral is regularized by the three- 
dimensional momentum cutoff A. Following Refs. j25il26ll . 
we introduce a parameter c as G\ = (1 — c)G+ and G2 = 
cG+, where < c < 0.5 and G+ > 0. The present model 
thus has four parameters of too, A, G+ and c. Assuming 
rno = 5.5 MeV, we have determined A and G+ from the pion 
decay constant f„ = 93 MeV and the pion mass M % = 138 
MeV at vacuum. Although c is an unknown parameter, we set 
c = 0.2 here, since it is known from the model analysis on the 
r\ — if splitting that c w 0.2 is favorable l32tl . 

For finite 8, parity is broken explicitly, so it is not a good 
quantum number anymore. Hence P-even and P-odd modes 
are mixed with each other for each meson. The "pion" mass 
M,,. is defined by the lowest pole mass of the inverse prop- 
agator in the isovector channel. It agrees with the ordinary 
pion mass when 8 = 0. Under the random phase approxima- 
tion lO . the inverse propagator is described by 



with 



det[l - 2GII{M 2 )} = 0, 



where 



G_ 
G+ 

n ss ( q 2 ) n SP { q 2 ) 
n PS ( q 2 ) n pp ( q 2 ) 



G = 



H{q 2 ) — ( ttPS(„i.. rri'i 



(28) 
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with 



n pp = AN f N c h-2N f N c {q 2 -AN 2 )I 2 {q 2 ), (31) 
n ss = AN } N c h-2N f N c {q 2 -4M 2 )I 2 (q 2 ), (32) 
n SP = n PS = —&NfN c M NI 2 (q 2 ) , (33) 
d A p 1 



h = 



(2tt) 4 p 2 - M 2 - N 2 ' 
4 - 1 



[(p±q/2) 2 - M 2 - N 2 ' 



(34) 
■(35) 



In this form, we can set a' = %' = 0, since we do not consider 
the isospin chemical potential. 

Applying the saddle-point approximation to the path inte- 
gral in the partition function, one can get 



(R(8)) 



detH n 



detH 



t_ e -/3V(n~n rc! ) 



(36) 



where (3 = 1 /T, Q (i7 re f) is the thermodynamic potential at 
the mean-field level in the original (reference) theory and H 
(-ff re f) is the Hessian matrix in the original (reference) theory 
defined by SB 



Hi 



d 2 n 



(37) 



For later convenience, the average reweighting factor (R(8)) 
is divided into two factors Ra and Rb- 



(R(6)) = RaRb 



(38) 



Ra = 



detP^ rc f 



detP 

Rb = e-W- *). 



(39) 
(40) 



For the TVj x N T lattice, the four-dimensional volume f3V is 
obtained by 



(41) 



Here we consider N X /N T = 4 as a typical example, following 

Refs. mm. 

We consider the following reference theory that has no sign 
problem: 

C = q'M + m rcf (8))q' - G + [{q'q 1 ) 2 + {q'l^rq') 2 ] 
- G_ [(q'Tq'f + (q'l l5 q') 2 } . (42) 

Here m lc f (8) is #-even mass defined below. We consider three 
examples as m Te {(9). 



A. The first example is reference A defined by 

m re f(8) = toa(6*) = to cos(6'/2). 



(43) 



In this case, the P-odd mass is simply neglected from 
the original Lagrangian ( TP7I ). 

B. The second example is reference B defined by 

m rc[ (8) = tob(0) 

= to cos(6I/2) + - {m sin(6»/2)} 2 . (44) 
a 

In this case, we have added the to-o -order correction 
due to the P-odd quark mass. Here a is a parameter 
with mass dimension, so we simply choose a = M„. 
The coefficient of the correction term is m^/M^ = 
0.129 MeV. 

C. The third case is reference C defined by 

TOref(0) = m c {8) 

= mo cos(8/2) + ^^- S m 2 (8/2), (45) 
v' 

This case also has the TOQ-order correction, but a is dif- 
ferent from reference B. The coefficient of the correc- 
tion term is m M 2 /M 2 , =0.114 MeV. 

Reference C is justified as follows. The pion mass M n (8) 
at finite 8 is estimated from the chiral Lagrangian as j23ll : 

M 2 (8) = !^|cos(0/2)| + ^f^sin 2 (0/2),(46) 



f2 

■I 7T 



J n 



where cto is the chiral condensate at T = 8 = and the coef- 
ficient I7 is evaluated by the 1/7V C expansion as 
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Fig. 1: 6 dependence of (a) the average reweighting factor and (b) 
M n at T = 100 MeV for the case of reference A. 



Fig. 3: 9 dependence of (a) the average reweighting factor and (b) 
M n at T = 100 MeV for the case of reference C. 
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Fig. 2: 6 dependence of (a) the average reweighting factor and (b) 
M n at T = 100 MeV for the case of reference B. 



The right-hand side of d46b is reduced to 

m M% 



l£o| 
f 2 



mn|cos(0/2)| 



Ml, 



\e/2) 



(48) 



Equation d48l i supports (|46*] i. 



III. NUMERICAL RESULTS 

If some reference system satisfies the condition that 
(R(9)) w 1, one can say that the reference system is good. 
As a typical example of the condition, we consider the case 
of 0.5 < (R(0)) < 2. This condition seems to be the mini- 
mum requirement. The discussion made below is not changed 
qualitatively, even if one takes a stronger condition. 

First we consider reference A. Figure 02 a) shows 8 depen- 
dence of (R(9)) at T = 100 MeV. The solid line stands 
for (R(6)), while the dashed (dotted) line corresponds to Ra 
(i?e)- This temperature is lower than the chiral transition tem- 
perature in the original theory that is 206 MeV at 9 = and 
194 MeV at 6 = n. As 9 increases from zero, (R(9)) also 
increases and exceeds 2 at 9 « 1.2. Reference A is thus good 
for 9 < 1.2. The increase of (R(9)) stems from that of i?e 
that depends on T. This means that the reliable 9 region in 
which 0.5 < (R(9)) < 2 becomes large as T increases. 

Figure Qlb) shows 9 dependence of M w at T = 100 MeV. 
The solid (dashed) line denotes M n in the original (reference 



A) system. At 9 = w, M T is finite in the original system, but 
zero in reference A. As a consequence of this property, Ra 
and (R(9)) vanish at 9 = it; see Fig. [Tf a )- This indicates that 
reference A breaks down at 9 = ir, independently of T. 

The same analysis is made for reference B in Fig. [2] As 
shown in panel (b), M n in reference B well reproduces that in 
the original theory for any 9. As shown in panel (a), however, 
the reliable 9 region in which 0.5 < (R(9)) < 2 is located 
only at 9 < 1.3. Therefore reference B is still not good. 

Finally we consider reference C. As shown in Fig. [3}b), 
in reference C well simulates that of the original theory 
at small and intermediate 9, but the former slightly underesti- 
mates the latter at large 9 around ir. As shown in Fig. 0b), 
however, (R(9)) satisfies the condition 0.5 < (R{9)) < 2 for 
all 9. Therefore we can think that reference C is a good refer- 
ence system for any 9. This is true for any temperature larger 
than 100 MeV. 



IV. SUMMARY 

We have proposed a practical way of circumventing the sign 
problem in LQCD simulations with finite 9. This method is 
the reweighting method for the transformed Lagrangian (0. 
In the Lagrangian, The sign problem is minimized, since the 
P-odd mass is much smaller than the dynamical quark mass 
of order tIqcd- Another key is to find out which kind of ref- 
erence system satisfies the condition (R(9)) 1. For this 
purpose, we have estimated (R(9)) by using the two-flavor 
NJL model and eventually found that reference C is a good 
reference system in the reweighting method. 
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